function CrackSIF
global glb

if all(glb.crack.state==0 | glb.crack.state==-1)
    glb.crack.KI(1,:)=0;
    return;
end

switch glb.solver.sif;
    case {1,'ansys'} % ansys
        CallAnsys;
    case {2,'abaqus'} % abaqus
        CallAbaqus
    case 3 % calculate
        CallMatlab
end
